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the technique as insensitive as possible to unknown or fluctuating 
phase and amplitude changes in the channel transfer function 
between the transmitter and receiver, frequency domain differential 
encoding techniques have been developed and their error performance 
calculated. It is shown that in the range of 2 to 8 bits/Hz of 
channel bandwidth efficiencies, differential encoding results in a 
penalty of 3 to 5 db in required E b /N 0 when compared to fully 
coherent multi-frequency modulation. Design procedures are 
presented that provide near optimum QAM constellations for fully 
differential coding and for the hybrid scheme of differential phase 
coding and absolute coding of amplitudes. 
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I INTRODUCTION 



Multi-frequency Modulation 

Multi-frequency modulation (MFM) is a digital modulation 
technique that is being researched at the Naval Postgraduate School 
for application in military digital communications systems. It is 
a high-da ta-rate, highly-bandwidth-ef f icient , modulation method 
with important potential for point-to-point, broadcast and packet 
switched network digital communications links. The method makes 
extensive use of digital signal processing techniques for its 
implementation . 

MFM is also known as orthogonal frequency division 
multiplexing (OFDM) , because of the orthogonality of the 
multiplicity of carrier frequencies that are used to transport the 
digital data. OFDM has been discussed in the literature in the 
context of data modems for the public switched telephone network 
(PSTN) , as a groupband modem, as a mobile HF modem, and for 
satellite broadcast of digital audio programs (DAB) [ 1 ] , [ 2 ] , [ 3 ] , [ 4 ] . 
One of the difficulties in the practical application of MFM (OFDM) , 
is that of carrier recovery for demodulation. Unlike ordinary 
bandpass digital modulation techniques that employ a single 
carrier, MFM uses a great number, perhaps hundreds or thousands, of 
carriers simultaneously. If the channel contains any envelope 
delay, that is phase shift that is not linear with frequency, then 
the channel transfer function must be known, or estimated, 
accurately if the carrier phases are to be used to carry data. 



During the past year we have concentrated on solving this 
problem. A solution that we have found works well in a practical 
channel [5], is to differentially encode information in the phase 
of adjacent tones in the same baud. The purpose of this report is 
to document the theoretical performance of this and of more general 
classes of MFM differential modulation. 

The underlying principles of MFM generation and demodulation 
were reported last year [6], The results of that report serve as a 
point of departure for this work. An MFM packet is shown in Figure 




1. The terms used in this report are consistent with those of [6]. 

Signal-to-Noise Ratio and E b /N 0 

It was shown in our previous report [6, pg 24] that for any 
given baud corrupted by additive white gaussian noise (AWGN) , the 
output of the DFT in the MFM receiver is the complex frequency 
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domain sequence 

Y(k) = S(k) + W(k) ; k 1 < k < k 2 (1.1) 

where the received signal baud is given by 



S(k) = A k cos0 k +jA k sin0 k , 



( 1 . 2 ) 



and 

A k =k x (2P k )*/2 

with P k the received tone power and k x is the number of points in 

the DFT. W(k) is a white zero mean complex gaussian frequency 

domain sequence with uncorrelated real and imaginary parts that 

2 

have equal variances k x AfN 0 /4. The frequency spacing between tones 
is if. By dividing (1.1) by the standard deviation of the noise 
components one obtains the unit variance complex sequence 



F(k) = R(k) +jl (k) = S'(k) + W'(k) 



(1.3) 



with mean value 



E[F (k) ] = S'(k) = A' k cos0 k +jA' k sin0 k (1.4) 

where 

A' k «(2E k /N 0 )* . (1.5) 

We have made use of the fact that P k /Af=E k is the received energy 
in the kth tone. This is because Af is the reciprocal of the baud 
length AT. 
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In some applications, the channel signal-to-noise ratio is the 
parameter most pertinent to system performance. We shall refer to 
the signal-to-noise ratio in the signal band as the wideband 
signal-to-noise ratio, SNRy B , and define it as the average signal 
power in the band to the average noise power in the same band. Thus 

SN^ = KP/N 0 W = E/N 0 (1.6) 

where P and E are the average tone power and tone energy 
respectively, and K is the number of tones in the signal band of 
width W =KAf. Furthermore, if there are Q bits transmitted per tone 
then E b = E/Q is the average energy per bit transmitted and 

Eb/N 0 = SNRy B / Q = E/QN 0 . (1.7) 

We shall use (1.5) through (1.7) freely when making performance 
comparisons of various types of MFM modulation. 

MFM Constellations 

Three classes of modulation will be considered in this report; 
phase modulation (PSK) , amplitude modulation (ASK) and combined 
amplitude and phase modulation. We shall refer to the latter, 
hopefully unambiguously, as QAM. (Ambiguity could arise from the 
fact the the PSK is coded and decoded using its quadrature 
components, but the context of our discussions should make the 
distinction between QAM and straight PSK clear.) In each case we 
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shall present results for both differential and non-differential 
encoding. A particuarily useful digital modulation format is a 
hybrid: differential phase coding and non-differential amplitude 
coding. The hybrid constellation for MFM has been developed by 
Mercury Digital Communications, Inc. [7], 

Differential modulation is obtained by encoding the input data 
symbols as changes in signal state occupancy of adjacent tone 
constellations in such a way that the changes are invariant under 
a rotation of the constellations. The principal purpose for 
differential encoding is to render the transmission insensitive to 
uncompensated channel phase shift. In MFM the tones can be placed 
very close together. If adjacent tones are close enough together 
then their uncompensated phase, that is the rotation of their 
constellations, will be nearly identical, and it cancels out in the 
decoding. 
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II ERROR PATTERNS AND PERFORMANCE MEASURES 



Symbols and Bits 

In MFM, Q bit symbols are coded into the phase and amplitude 
of each of the multiplicity of carrier frequency tones during each 
baud. The symbols are partitioned into M bits that are coded into 
tone phase and N bits that are coded into tone amplitude (Q=M+N) . 
We adopt the notation 

s = p,a (2.1) 

to designate the symbol and its partition into phase and amplitude 
bits. For example, with 4 bits per tone, one may have 



1010 = 101,0 ( 2 . 2 ) 

and three bits of phase, in this example 101, are coded as one of 
eight possible tone phases and one bit of amplitude, in this 
example 0, is coded as one of two possible tone amplitudes. 

Errors may occur in the receiver in any of the Q bit 
positions. An error pattern is defined as a Q bit symbol 

e = ep,e a (2.3) 

which we also partition into phase bit errors and amplitude bit 

Q M 

errors. There are 2 total error patterns made up from 2 phase 
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error patterns for each of the 2 amplitude error patterns. 

The output of the receiver is the decoded symbol 

d = s + e (2.4) 

where the plus sign designates bit by bit binary addition (logical 
"or") of the transmit symbol with the error pattern. The decoded 
symbol is also partitioned into decoded phase and amplitude bits, 

d = d p ,d a (2.5) 

A symbol error occurs anytime that d * s, that is for any of the 
2 Q -1 non-zero error patterns, e * 0. 



Symbol Error Probabilities and Bit Error Rate CBER) 

The probability of symbol error is 

P s = S Pr [ e*0 | s=s ; ] Pr [ s=s ; ] . (2.6) 

« 

If we assume proper source coding then each of the 2 ° data symbols 
will be equally probable and 

P s = 2 ° S Pr[e*0 | s=Sj] . (2.7) 

t 

Performance of digital communications systems is customarily 
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specified is terras of bit error rate (BER) . BER is defined as the 
expected number of bit errors in a very long sequence of symbols 
divided by the number of bits in the sequence. For example, if the 
expected number of bit errors in a sequence of one million bits is 
one, then the BER is 10 6 . As long as symbols that make up a 
sequence of bits are statistically independent and the error 
probabilities are stationary, then the BER is also the expected 
number of bit errors in any given received symbol divided by Q, the 
number of bits in the symbol. We shall adopt this definition of BER 
in this report. 

Let B be the number of symbol bits in error, that is the 
number of ones in e. Then 



E[B] = £ E[B| s=£,]Pr [£=£,] 



( 2 . 8 ) 



and again assuming equally probable symbols, 

(2.9) 

( 2 . 10 ) 



-o 



E[B] =2 £ E [ B | s=s j ] 



with 



BER = E[B)/Q 



in either case. 

Single and Multiple Bit Errors 

BER is a normalized expected value. Sometimes it is most 
convienent to compute it as a weighted superposition of terms 
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corresponding to single bit errors, double bit errors, etc, up to 
Q bit errors. That is, 

q <5 

E [ B] = 2 b Pr [B=b] =2 2 {2 b Pr [ B=b | s=s- ] } (2.11) 

\ v \ 

where, Pr[B=b|s=Sj] are the conditional bit error probabilities. 
These can be computed from the conditional error pattern 
probabilities as follows; 

Pr [ B=b | s=s,- ] = 2 Pr[e|e contains b ones,s=Sj] (2.12) 

where the sum is over the N b = (Q,b) error patterns containing b 
ones out of Q bits. For example, there are N 1 = (Q,l) = Q single 
bit error patterns, N 2 = (Q,2) = Q(Q-l)/2 double bit error patterns 
etc. 

In summary, the symbol error probability is computed from 
(2.7) and the BER from (2 . 10) , (2 . 11) and (2.12). Both computations 
require the conditional error pattern probabilities Pr [ e=e - 1 s=s ■ ] . 
The objective of the next three Chapters is to obtain theoretical 
expressions for these for various MFM and MFDM constellations. 
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Ill PHASE MODULATION 



In this Chapter MFM and MFDM systems that use constant 
amplitudes and employ only phase modulation to encode the Q symbol 
bits (M=Q, N=0) will be considered. Two types of modulation will be 
analyzed; coherent modulation/demodulation and differential 
modulation/demodulation. In coherent modulation, the absolute phase 
of each of the carrier tones must be known at the receiver. This 
form of modulation should only be used in a linear phase channel, 
or in a channel that has been equalized to linear phase. Any 
residual phase jitter will affect performance as noise, and any 
significant unequalized non-linear phase (group delay) will most 
probably render coherent phase modualtion unworkable at any SNR. 
However, its analysis, assuming perfect phase equalization, is 
included for comparison with differential modulation. 

Differential phase modultion of MFM, MFDM, encodes phase bits 
as the phase difference between adjacent carrier frequency tones in 
the same baud. Although it is not necessary that adjacent tones be 
used, or that the same baud be used (we could encode using the 
phase difference of each of the same tones on successive bauds as 
is done in conventional single carrier differential phase 
modulation (DPSK) transmission) , the theory behind differential 
encoding in the frequency domain using adjacent tones is that the 
differential phase jitter between adjacent tones introduced by the 
channel can be reduced to any desired level by putting the tones 
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sufficiently close together in frequency. Recall that the tone 
spacing in MFM is Af = 1/AT. Therefore differential phase shift 
through an arbitrary unequalized channel is inversely proportional 
to baud length. 

Coherent Phase Modulation/Demodulation 

A phase constellation for 2 M -ary PSK is generated by utilizing 
a constant tone amplitude and 2 phases. The phases are equally 
spaced at intervals A0= 2ir/2 H , then rotated by adding A0/2. Thus 
the modulation phase bits are encoded according to the 
transformation, 



Elkl i - d0(k) = i(k)A0 



<p (k) = d0(k) + A0/2 



i = 0,1, 2-1 



(3.1) 



Such a constellation for 8 phases is shown in Figure 2. 

The maximum likelihood demodulation rule (MLD) is to decode 
the complex point Y(k) generated by the receiver DFT(the matched 
filter receiver for MFM, see [6]) to the nearest point of the 
constellation. The decoding boundaries for MLD are sectors of 
anglular width A <p as illustrated for an 2 M of eight in Figure 2. We 
shall refer to this type of decoding as "sector decoding" (SD) . 
Notice that SD is independent of tone amplitude A, or equivalently 
of E/N 0 . This makes the implementation of SD in the receiver 
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independent of system gain, a highly desireable feature. Of course 
system performance is not independent of SNR. 



L 




Figure 2 An Eight PSK Constellation with Sector Decoding 
Boundaries 



a. Symbol Error Probability 

Phase symbol error decoding probabilities are the same for 

M 

each of the 2 symbols. This is so because the statistical 
properties of the complex noise W(k) are invariant to phase 
rotations. That is, the real and imaginary parts of W(k)exp(j0) 
are, like those of W(k) , uncorrelated, equal variance gaussian 
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random variables. Consequently, the probability of symbol decoding 
error given by (2.7) is the same as the conditional probability of 
decoding error 



P s = Pr [e*0 1 s=Sj ] , for any i. (3.2) 

The probability that e # 0. when s = s,- is the probability that 
Y(k), or equivalently F(k), is not in the sector of angular width 
A 0 defined by the angles iA0 and (i+l)A0. 

i.) QPSK (M = 2) 

We include this case for completeness as the results are 
well known and were presented in the previous report [6]. Also, for 
more phase states and narrower phase sectors, the results that will 
be presented in the next section are only approximate (the 
approximation becomes more accurate as the sectors become more 
narrow) . Exact results are easily found for QPSK as follows: 

The probability that F(k) is not in the first quadrant, 
given ^ is simply one minus the probability that F(k) is in the 
first quadrant. The probability that F(k) is in the first quadrant 
is the probability that R(k)>0 and that l(k)>0. Since the real and 
imaginary parts of F(k) are statistically independent; 

P s = 1 - Pr{R(k) >0}Pr{I (k) >0} . (3.3) 
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These probabilities are easily found using the fact that R(k) and 
I(k) are unit variance gaussian random variables with mean values 
determined from (1.4) -(1.6) to be 

E{R(k)> = E ( I (k) } = A' k /(2)* = P (3.4) 

with 

p = { SNRy B ) % ={E/N 0 }* , 

where we have used the fact that 0 k = -n / 4 for s = Sg, the symbol 
assumed as transmitted on tone k. Consequently, 

P s = 1 -[ 1 - Q(P) ] 2 * 2Q(p) (3.5) 

is the symbol error probability. The approximation is very good for 
SNRyg greater than about 10 db. Note from (1.7) that for QPSK 

p 2 = E/N 0 = SNRyg = 2E b /N 0 . (3.6) 



ii.) 2 M -ary PSK 

Let p(A9) be the probability that F(k) falls outside 
an angular sector of width ±A©/2 with respect to <p s , the angle of 
F(k) given the transmission of s,- in the absence of noise. That is, 

p(A0) = Pr[|a - 0j | >A0/2 ] (3.7) 
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where a is the angle of Y(k) (or equivalently of F(k)). 

It is known that for A© < m/2, this probability is well 

approximated by ([ 8 ]; Hakin pg 317). 

p(Ae) * 2Q( A ' k sin ( A 0/2 ) } = 2Q{ (2 ) % psin (A 0/2 ) } (3.8) 

and that the approximation becomes increasingly more accurate with 
increasing p and decreasing A©. A PSK symbol error occurs if 
| a - (p ] \>k<p/2. Therefore, the PSK symbol error probability is 

P s = p(A0) * 2Q{ ( 2 )*psin(A 0 / 2 ) ) . (3.9) 

We note that (3.9) with A <p = m/2 , which is QPSK, agrees with the 
approximation generated by binomial expansion of the exact equation 
derived for QPSK and presented in (3.5). 



Differential Phase Modulation/Demodulation 

In MFDPSK, the phase bits for tone k are encoded into equally 
spaced phases 

p(k) ,• -*• d 0 (k) = i(k)A 0 ; i=l, 2 , . . . . 2 M - 1 , (3.10) 

and this differential phase shift is added to the phase of the 
previous tone to obtain the tone phase of tone k as 



15 



0(k) = 0(k-l) + d0(k) . 



( 3 . 11 ) 



An initial condition must be chosen. By choosing the phase of the 
first carrier frequency tone as 

0(k.,) = A0/ 2 (3.12) 

the constellations of all tones are rotated by A0/2 thereby 
yeilding transmit tone phases identical to those used in coherent 
MFPSK discussed in the previous section. 

Differential decoding the phase bits is accomplished by 
computing the complex rotated product 

D (k) = F (k) F* (k-1) exp[ j A0/2 ] ; k.,+l<k<k 2 (3.13) 

for adjacent digital frequencies in each MFDM baud. In the absence 
of noise, the phase of D(k) is the phase difference between tone k 
and tone k-1 plus A0/2, that is i(k)A0 + A0/2. Comparison with 
(3.1) reveals that the complex rotated products have exactly the 
same constellations as described for coherent PSK. The differential 
receiver optimally decodes the phase bits at this stage by using 
the same sector decoding rules described above, 
a. Statistical properties of the Complex Rotated Product. 

Define the complex random product 



F p (k) = F (k) F* (k-1) = S' (k) S ' * (k-1) + W p (k) (3.14) 
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with 



w p (k) = W' (k)W'*(k-l) + S' (k)W'*(k-l) + W' (k)S'*(k-l) . (3.15) 



Since W' (k) is zero mean, it follows directly that W p (k) has zero 
mean and therefore that the mean of F p (k) is 

E[F p (k)] = S' (k)S'*(k-l) = {2(E k E k . 1 )VN o }exp[jd0(k) ]. (3.16) 

The real and imaginary parts of the complex product are 

R p (k) = Re [ F (k) F* (k-1) ] = R(k)R(k-l) + I(k)I(k-l) 

(3.17) 

I p (k) = Im[F(k) F*(k-1) ] = I (k) R(k-l) - R(k)I(k-l). 

Because of the statistical independence between adjacent tones and 
between the real and imaginary parts of the same tone, it follows 
directly that 

Var{Re[W p (k) ] }= Var{R p (k) } = Var [R(k) R(k-l) ] + Var [ I (k) I (k-1) ] 

(3 . 18) 

Var { Im[W p (k) ] }= Var{I p (k)} = Var [I (k)R(k-l) ] + Var [R(k) I (k-1) ] . 

Let X and Y be statistically independent, non-zero mean random 
variables. Then 
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Var [XY] = E[X 2 Y 2 ] - {E[XY]} 2 = 

E[X 2 ]E[Y 2 ] - { E [ X 3 } 2 { E [ Y ] } 2 = 

a x 2 a y 2 +{E[X] } 2 a y 2 + {E[Y]} 2 a x 2 . (3.19) 

It can be seen that the variance of a product contains a noise 
times noise term and two signal times noise terms; the later two 
are due to the mean values of X and Y. Recalling now that the real 
and imaginary values R(k) and I(k) have been normalized to unit 
variance for all k, then it follows that 

o p 2 = Var{R p (k) } = 2(1 + E k /N 0 + E k . 1 /N 0 ) 

(3.20) 

a p 2 = Var{I p (k) } = 2(1 + E k /N 0 + E k _ 1 /N 0 ) . 

Here, the first term is due to noise times noise, the second two 
are the signal times noise terms. This is important because at high 
SNR, the signal times noise terms dominate and R p amd I p are 
approximately gaussian with mean values given by (3.16) and 
identical variances given by (3.20). 

It is important to determine if the real and imaginary parts 
of the complex product noise terms, W p (k) are correlated. We note 
that, because W p (k) is zero mean, 

Covar{Rp(k)I p (k) ) = E{Re[W p (k) ] Im[W p (k) ] } . (3.21) 
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The covariance is defined as 



Covar{Rp(k)I p (k) } = E{R p (k) I p (k) } - E{Rp(k) }E{ I p (k) } . (3.22) 

which we find by using (3.17) and the previously cited properties 
of the noise to be identically zero (this is so at all SNRs and is 
not dependent on the gaussian approximation) . Zero covariance 
insures that D(k) , the rotated complex product of (3.13) also has 
uncorrelated real and imaginary parts with equal variances as given 
by (3.20). At high SNRs it tends to the gaussian. 

Dividing (3.13) by a p produces a normalized differentially 
demodulated complex random variable D'(k) with unit variance real 
and imaginary parts and mean value given by, 

E [D • (k) ] = (S' (k)S'*(k-l)/a p }exp( jA0/2) = 

M' k cos(d0(k)+A0/2) +jM' k sin(d0(k)+A0/2) (3.23) 

where 



M 'k =[( 2E k E k-l)VN 0 ]/[l+E k /N 0 +E k . 1 /N 0 ) ,/2 . (3.24) 

Comparing (3.23) and (3.24), with (1.4) and (1.5) it can be seen 
that differential modulation causes a reduction in SNR for 
decoding. If, for example, adjacent tones have equal energy and if 
the SNRyg = E k /N 0 is much larger than one, then 
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(3.25) 



M\ * (E k /N 0 } % =p k 

which is 3 db less than (1.5) predicts for the coherent PSK. 

The symbol error probabilites are determined as before. So 
long as the SNR is sufficiently high to use the gaussian 
approximation, then the probability of a symbol error for MFDQPSK 
is 



P s = 1 - [1 - Q{M' k /(2)*}] 2 ~ 2Q{M' k /(2)V (3.26) 

The probability that D(k) falls outside a sector of width A© 
is 



p (A©) ~ 2Q{M' k sin(A©/2) } (3.27) 

so the probability of a differential symbol decoding error in 2 - 
ary MFPSK is 

P s = p(A <p) (3.28) 

with M' k given by (3.24) and A 0 = 27T/2 M . 

Bit Error Rates for MFPSK and MFDPSK 

In order to compute BERs, the expected number of bit errors 
per symbol must be computed and divided by the number of bits per 
symbol (see (2.10)). The probability of single, double, triple bit 
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errors etc. depends on the transformation that is used , that is it 
depends on (3.1), (3.10), 



Eikli 



i (k) A0. 



(3.29) 



Because the most probable errors are in the sectors that are 
adjacent to the sector containing the transmitted symbol, and the 
next most probable are in the next set of adjacent sectors, etc., 
we have elected to use a reflex or Grey code to code the phase 
bits. In this way, adjacent sectors only differ by one bit, so 
that the most probable symbol decoding errors are single bit 
errors, next most adjacent sectors differ by two bits, so they 
yield double bit errors. In this way, we hope to minimize BER when 
the errors are due to AWGN. 



The 


Grey reflex code for 2, 


3 , and 


4 phase bits 


is given 


reference 


in Table 


3,1 










P.IK).,- 






i(k) 




M=2 


M=3 


M=4 


M=2 


M=3 


M=4 


00 


00,0 


00,00 


0 


0 


0 






00,01 






1 




00,1 


00,11 




1 


2 






00,10 






3 


01 


01,1 


01,10 


1 


2 


4 






01,11 






5 
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01,0 


01,01 




3 


6 






01,00 






7 


11 


11,0 


11,00 


2 


4 


8 






11,01 






9 




11,1 


11,11 




5 


10 






11,10 






11 


10 


10,1 


10,10 


3 


6 


12 






10,11 






13 




10,0 


10,01 




7 


14 






10,00 
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TABLE 3.1 Reflex Coded Phase Bits 

Note from Figure 2 that the two most significant bits are 
determined soley from the quadrant in which the received complex 
point occurs. The remaining bits depend only on the sector occupied 
by the absolute values of the real and imaginary parts of the 
complex values. These observations are useful for designing 
encoding and decoding algorithims for variable numbers of bits of 
phase modulation in a systematic fashion (see, Mercury Digital 
Communications patent application [7]). 

a. MFQPSK and MFDQPSK Bit Error Rates. 

It is easy to see that the two quad bits are decoded as 

follows ; 
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Most Significant Bit -* 0 if I(k) > 0; 
Most Significant Bit -*• 1 if I(k) < 0. 



Next Most Significant Bit -*■ 0 if R(k) > 0; 
Next Most Significant Bit -+ 1 if R(k) < 0. 



Here R(k) and I(k) are the real and imaginary parts of Y(k) or of 
D(k) for coherent MFQPSK or MFDQPSK respectively. Because, in both 
cases they are statistically independent random variables, then 
QPSK is equivalent to sending two independent one bit symbols and 
consequently the BER is just the symbol (bit) error probability of 
either one. Thus, for QPSK 



BER = Q (x) 



(3.30) 



where 



for Coherent MFQPSK; 



x = A' k /(2)* = P 



(3.31) 



and for MFDQPSK; 



x = M' k /(2) 



= P/2 



(3.32) 
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b. 2 M -ary Coherent MFPSK and MFDPSK. 

Due to the symmetry of the phase sectors and the Grey 
coding algorithim employed in the PSK and DPSK, it is clear that 
the conditional probabilities of B bit errors in a symbol of M 
phase bits are all equal. That is, Pr[B=b|s =s,] is the same for 
all s,. . and therefore so is E[B|s = s,]. Thus, for PSK systems 

BER = E[B|s]/M (3.33) 



and 



E[B|s] = l(Pr[B=l|s]] + 2 { Pr [ B=2 | s ] } + ...+ M( Pr [ B=M | s ] } . (3.34) 

Figure 3 illustrates B , the number of bit errors for M =4 in 
each of the 16 sectors with respect to the transmitted symbol s = 
00,00. Note that of the 4 single bit errors, 2 are in adjacent 
sectors, and of the 6 double bit errors, 2 are in sectors adjacent 
to the single bit sectors. We shall assume for the remainder of 
this report that the probability a received complex value is more 
than two sectors away from the transmit sector is negligible 
compared with the probability of being within the two closest 
sectors. Thus, 



Pr [B=l | s] a p(A0) - p(3 A0) 



(3.35) 
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Pr [B=2 | s] ~ p(3A0) 



and consequently 



BER *{p(A<*>) + p(3A0)}/M. 



(3.36) 



with p(A0) determined by (3.8) for coherent MFPSK and by (3.27) for 
MFDQPSK. 



Equations (3.30) 
and (3.36) represent the 
principal results of 
Chapter HI. They are 
used to compute BER for 
coherent phase coded and 
for differentially phase 
coded MFM. Grey coding 
of the phase bits into 
the phase sectors has 
been assumed. 




Figure 3 Distribution of Bit Errors B for 
s =00 , 00. 
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IV AMPLITUDE MODULATION 



In this Chapter, results are presented for digitally coding 
information into the amplitude of MFM carrier tones. For the non- 
differential case, these are but a recitation of well known results 
for incoherent amplitude shift keying since all the tones are 
orthogonal and are in effect processed in parallel. The advantage 
of amplitude modulation is that amplitude coding does not require 
any knowledge of the individual carrier tone phases for 
demodulation. It does require knowledge or control of the gain of 
the channel at each of the carrier tones in order to set the 
thresholds for the decision regions. Consequently, a technique for 
differentially encoding data as the change in amplitude between 
adjacent carrier tones will also be described and analyzed. This 
method and an apparatus for its implementation have been described 
in a US Patent Application [7 ] . So long as the channel gain is 
nearly identical at adjacent frequencies, the differential 
amplitude coding system is independent of both the phase and the 
amplitude response of the channel. 

Multilevel Amplitude Modulation 

Multilevel, multi-frequency amplitude modulation encodes N 
amplitude bits into one of 2 N amplitude levels in each of the 
carrier tones. The coding is prescribed by a transformation; 

ajklj - A k = i(k)AA + A^ ; i = 0,1...2 N -1, (4.1) 
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and as with the MFPSK coding, a Grey code is used so that adjacent 
amplitude values correspond to only one bit difference in 
codewords . 

Given the transmission of codeword a (k) . the normalized 
decoded complex DFT coefficient F(k) is a unit variance gaussian 
complex random variable with uncorrelated real and imaginary parts 
and average values given by (1.5). The magnitude of F(k) , is a 
Rice-Nakagami random variable, R, with probability density function 



f R (r[ a(k) )= rexp{-(r 2 +A' k2 )/2)I 0 (rA' k ) (4.2) 

Now, for A' k »1, integrals of f R (r) in the vicinity of A' k can be 
approximated by integrals under the gaussian distribution with mean 
value A ' k and unit variance, (see e.g. [9], pg 504 and pg 625.) 

Amplitude bits will be decoded in accordance with the MLD rule 
as illustrated in Figure 4. Conditional amplitude symbol error 
probabilities are; 

PJA'k - * Q(4A'/2) 

M*'.ln <A 'k < » 2Q(4A'/2) (4.3) 

P s [ A ’k “ A ' nraxJ » Q(4A'/2) 

Assuming equally likely symbols, the average symbol errror 
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probability is 



P s = 2 N S P s [A' k ] =2(1- 2 N ) Q(AA'/2) (4.4) 




Figure 4 Multi-Level Amplitude Decoding 



Average SNR^.E/Np and E b /N 0 

When amplitude modulation is used, the carrier tones have 
different power. However, we would like to evaluate the error 
probabilites in terms of an average SNRy B and average 
Let A 1 0 =E[A k '] be the average signal amplitude of tone k and let 
the increment between amplitudes be proportional to the average 
amplitude, that is let AA'=aA' 0 . Then it is easily shown that the 
average tone energy E is given by 

E =E[E k ] = E 0 { 1 +(a/2) 2 (2 N -l) [l+(2 N -2)/3]} (4.5) 
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where 



E 0 = (A' 0 ) 2 N 0 /2 (4.6) 

is the energy associated with the average amplitude of tone k. We 
shall be most interested in zero, one, two and three bits of 
amplitude coding for which N is zero, one, two and three 
respectively and (4.5) shows that the average energy is related to 
the energy of the average amplitude tone by 



N=0 : E = E 0 

N=l : E = E 0 { 1 + (a/2 ) 2 ) 

(4.7) 

N=2 : E = E 0 (l +5 (a/2 ) 2 ) 

N-3 : E = E 0 {1 +21 (a/2 ) 2 ) . 

The average SNRy B and average Et/N 0 are related to the average tone 
energy E as in (1.6) and (1.7), viz; 

SNRy B = E/N 0 = p (4.8) 



and 



V N 0 = E /QN 0 . 



(4.9) 
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In order to evaluate symbol error probabilites , we need to 
evaluate Q (AA '/2 ) ♦ Now recalling that AA' has been defined as a 
fraction of the average tone amplitude, then it follows that 

Q(AA'/2)=Q(aAV2) = Q (a { 2E C /N 0 } % /2 ) 

= Q(<2E/N 0 }*f(a,N)) (4.10) 

with 

f(a,N) = a/[2{l+(a/2) 2 (2 N -l) [l+(2 N -2)/3] }*] (4.11) 

It is clear that the function f(a,N) is of central importance 
for controlling error probabilites. Notice that there is a maximum 
value for a because the smallest amplitude level must be non- 
negative and that the maximum value decreases with increasing N. 
That is. 



0 < a < 2/ (2 N -1) . (4.12) 

We note that f(a,N) can be reasonably approximated by a/2 for 
values of a less than one. It is also clear that if one is just 
using amplitude modulation, a should be made as large as possible, 
which is the same as spreading the tone amplitude levels as much as 
possible, in order to minimize error probabilities. Consequently, 
with amplitude modulation alone, performance is limited by 

f (<W N > = (1/(2 N -1) }{1/[1+(1/(2 N -1)) (l+(2 N -2)/3) ]*} (4.13) 
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For one, two, three and four bits of amplitude modulation, (4.13) 
is given in Table 4.2 below. 



N 



f-tow^N 



1 

2 

3 



1/2 

v 

1/14 

1/70 



% 






4 1/310 

TABLE 4.2 

The Maximum of the Spread Function 



Finally, combining (4.4) and (4.10), the symbol error 
probability for amplitude modualtion alone of the tones is given by 

P s = 2(1-2' N )Q{ (2)V(a,N) ). (4.14) 

Amplitude and phase modulation symbol error probabilites for 
different numbers of bits may be compared by comparing (4.14) with 
(3.9) . ( It should be noted that for the situation of a = cc max , that 
the minimum amplitude A' m1n =0 and therefore the probability 
distribution is Rayliegh, not Gaussian as assumed in (4.3). Since 
we are primarily interested in combined amplitude and phase 
modulation constellations for MFM for which 0 amplitude tones are 
not permissible, we shall not be concerned with presenting the 
necessary correction here.) 
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BER's for Amplitude Modulation 



Bit error rates for the amplitude bits are determined as 
presented in Chapter II. That is. 



We shall approximate (4.16) using single and double bit 
errors. This is a good approximation by virtue of the Grey coding 
of the amplitudes bits. Thus, assuming equally likely amplitude 
codewords, 



These conditional error probabilities depend on the whether, the 
codeword a corresponds to A' mjn or A'^ or to an in-between value 
of A'. We are principally concerned with one, two and three bits 
of amplitude modulation. The evaluation of (4.16) for these cases 
is the following: 



BER = E[B]/N 



(4.15) 



with 



E[B] = 2b (2 Pr[B=b|a]Pr [a] } (4.16) 



E [ B) = 2 H (2 Pr [ B=1 | a] +22 Pr[B=2|a]} (4.16) 



N=l: 



E [ B] = Ps = Q(AA'/2) = BER 



(4.17) 



N=2 : 



E [ B] = 1.5Q(AA'/2) + Q(3AA'/2)= 2 BER 



(4.18) 
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N=3 : 



E[B] = (7/4)Q(AA'/2) + (3/2)Q(3AA'/2) = 3BER (4.19) 

with 

AA'/2 = { 2 } % pf (a , N) (4.20) 

Differential Amplitude Modulation 

In order to decode the amplitude bits of MFM carrier tones, it 
is necessary to establish accurate values for the channel gain at 
each carrier tone frequency. Although this is less difficult to do 
than to determine phase shift, it may still be of advantage to have 
a modulation scheme which does not require either phase or 
amplitude information in order to decode the symbols. In mobile 
systems with bandpass signals subject to fading this will almost 
certainly be the case. Also, in switched systems where a premium is 
on connect time, allocation of valuable connection time to 
measurement of channel gain may represent a limitation. 
Consequently, we have investigated the performance of MFM 
modulation with differentially encoded amplitude information. 
Although encoding more than one bit differentially is conceptually 
possible, to date we have only analyzed and tested modulation 
schemes involving one bit of differential amplitude coding where 
the coding is, as with the differential phase modulation, 
differential between adjacent tones. Thus, two amplitudes, A^ and 
A^ih that are separated by AA = aA 0 = Q:{A max + A min }/2 are used. 
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Differential encoding the amplitude bit of tone k is accomplished 
by making A k the same as A k .., if a=0 and changing A k to the amplitude 
different from the amplitude of A k _., if a=l. The decoding rule for 
a is illustrated by Figure 5. The crosshatched region is decoded as 
a=0, that is it is assumed that the amplitudes are the same in this 
region. The remaining region is decoded as a=l, that is it is 
assumed that the amplitudes differ in this region. This rule is 
described by the algorithim; 

"if {1/(1 -a/2) } |Y(k-l) | < | Y (k) |<( l-a/2 ) | Y (k-1) | then a=0 

else a=l" . 

Clearly the rule is 
independent of channel gain so 
long as adjacent tones have 
essentially identical gain. 

To find symbol and bit 
error probabilities, the 
conditional probability of |y| 
being outside the decision 
regions must be computed. Recall 
that we have established the 
mean and variance of |y| and the 
fact that adjacent magnitudes 
are statistically independent. Consider the dotted lines with 
slope of one passing thru the points (A 0 ,A mjn ) and (A mjn ,A 0 ) . We shall 




Decoding 



34 



use these as approximations for actual decision boundaries shown by 
the solid lines. Accordingly, it can be seen that the average 
symbol error probability which is also the BER is given by; 

BER = PS = (3/2)Q{ (l/2 % ) AA'/2} - ( 1/2 ) Q { ( 1/2 ) *3 AA ' /2 } . (4.21) 

Comparing (4.21) and (4.17), we note that a differentially encoded 
amplitude bit requires approximately 3db more signal-to-noise ratio 
than an absolutely encoded amplitude bit in order to obtain the 
same BER. 
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V. QUADRATURE AMPLITUDE MODULATION 



The purpose of this Chapter is to describe more complex 
modulation constellations, constellations that combine the 
amplitude and phase modulation schemes discussed in the previous 
two Chapters. As mentioned in the Introduction, we shall refer to 
these as "MFQAM systems" even though they will not have Z 2 (that is 
rectangular) grid points as is commonly found in single carrier QAM 
constellations . 

To be more specific, we will first consider systems that 
utilize differential phase modulation of M bits and one bit of 
differential amplitude modulation. These systems will be referred 
to as "fully differential systems". Also, we will consider systems 
that use differential phase coding and absolute amplitude coding. 
A system of this type with M bits of differential phase modulation 
and N bits of absolute amplitude modulation will be referred to as 
an "(M, N) Hybrid". 

Single Amplitude Bit MFQAM Systems 

These systems accept Q=M+1 bit codewords s which are 
partitioned into two quad bits that determine the phase difference 
of adjacent tones to the nearest ninety degrees , M-2 delphase bits 
that are coded into phase within the quadrant and one amplitude bit 
the is encoded as tone amplitude. Thus a codeword is divided into 
the three subsets 
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s = g,d,a . 



(5.1) 



As described in Chapters II and III, all coding is done in 
accordance with a Grey (reflex) code such that nearest neighbor 
points in QAM constellations differ in only one bit position. The 
BER for these constellations is determined from (2.10) using (2.11) 
and (2.12). We shall consider only single and double bit errors. 
Accordingly, 



are to be evaluated for all s,-. 
a. Single Bit Errors 

Single bit errors occur when an amplitude bit decoding error 
occurs and all the phase bits are decoded correctly or when a 
single bit phase error occurs and the amplitude bit is decoded 
correctly. 

These correspond to the error patterns; 



Pr[B=l|s=Sj] and Pr[B=2|s=s 1 ) 



(5.2) 



e = 0,1 and e = e p1 ,0 



(5.2) 



where e p1 is a single bit phase error pattern. 
Therefore the single bit error probability is. 



Pr [B=l] = Pr[e a =l,ep=0] + Prfe^e^O] (5.3) 
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If we assume amplitude and phase decoding errors are independent 
(which is certainly not strictly true but is probably a 
satisfactory approximation at low error rates) then, 

Pr [B=l] = Pr[e a =l]Pr[ep=0] + Pr [e,,, ] { 1-Pr [e a =l] } 

~ Pr[e a =l] + Prtfip,] (5.4) 

where the approximation is accurate at low symbol error rates. 

The probability of an amplitude bit error is clearly 

independent of the phase of tones and is given by (4.17) and (4.21) 
for absolute and differential amplitude coding respectively. 

The probability of phase bit errors is identical for all 

phases but depends on the amplitude of the complex rotated product 
which we have called M k and whose statistics are given in (3.24). 
A single amplitude bit, that is two amplitude levels, results in 
three different complex rotated product amplitude levels for M k 
depending on whether adjacent tones both have amplitudes A^, both 
have amplitudes A^, or one has amplitude A^ and the other 

amplitude A max . Since we have assumed all symbol patterns are 

equally likely, then the first two situations each occur with 
probability 1/4 and the latter with probability 1/2. 

Let us designate the three levels of M' k as M , (nin ,M , mid and M'^. 
Then the probability of any particular phase bit error pattern is 

Pr[e p ] = >iPr[e p |M'. id ] +(, ( Pr(e p |M'„, n + Pr^lM'^] ) (5.5) 



38 



The probability of single (and double) phase bit error patterns are 
given in (3.35) where p(A0) is given in terms of M' k by (3.27) for 
MFDPSK. The three levels for M' k are found from (3.24) for E/N 0 > 
1 to be, 



M'-in = <WN.)‘ 

H'.ld = (E.i„W EN .)‘ (5-6) 

In terms of the amplitude level spread parameter a, 

E m\n = E o k ( 1 ~ a / 2 )' E max^ = E 0 *( 1+a / 2 ) and E 0 ' the energy of an average 
amplitude tone is related to the average energy of the tones by 
(4.5). Specifically we find that, 

M' min = { (E/N 0 ) (l-a/2)/(l+a/2) }* 

M' mid = { (E/N 0 ) (l+(a/2) 2 ) )’ A (5.7) 

M 'max = {(E/N 0 )(l+a/2)/(l-a/2)}\ 

Thus (5.7) in conjuction with (5.5) and (3.35) determines the 
probability of single phase bit errors in terms of the p= SNRyg = 
E/N 0 , the number of phases 2 M and the parameter a. Finally the BER 
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due to single bit errors is well approximated by 



BER 1 =Q {Pr[e a =l] +Pr[e p1 ]} 



(5.8) 



with 



Q(pa 0 ) for hybrid coding 



Pr [e =1] = 



((3/2)Q(pa 0 /2*) -(l/2)Q(3pa 0 /2 % ) for fully 



differential coding, and 



Pr [Gpi ] = Q(pa 2 )-Q(3pa 2 )+^[Q(pa 1 )-Q(3pa 1 ) +Q ( pa 3 ) -Q (3pa 3 ) ] 

with 

p = (E/N 0 ) % 

a 0 = a/{2[l+(a/2) 2 ]}* 

a, = sin(7r/2 M ) { (l-a/2)/ (l+a/2) )' /! 

a 2 = sin(7r/2 M ) {l+(a/2) 2 ) % 



40 



a 3 = sin(7T/2 M ) { (l+a/2)/(l-a/2) >*' 



Examination of (5.8) shows that all the a values increase with a 
except a 1 which determines the phase errors due to complex products 
of the minimum amplitude tones. This makes sense as we know that 
the amplitude errors decrease with increasing amplitude difference 
and the large amplitude tones will have more energy with increasing 
amplitude differential. 

One design procedure we have investigated is to choose a such 
that a 1 = a 0 for the hybid, i.e. absolute amplitude bit coding and 
a, = a 0 /2 A for differential amplitude coding. This leads to the 
following values of a; 



Q 


M 


a Fullv Differential Codinq 


a Hybid 


Codinq a 


3 


2 


.938 


.728 


.88 


4 


3 


.593 


.447 


ID 

• 


5 


4 


. 335 


.250 


.32 


6 


5 


.180 


.130 


. 18 



TABLE 5.1 

Amplitude Spread Factor for 
a Single Bit Encoded as Amplitude 

From a purely geometrical point of view, values of 



41 



a=A0/ (1+A0/2) space the constellation points approximately 
equally. These values are listed for comparison in Table 5.1 too. 
Using the optimum values of a yeilds the a coefficients shown in 
Table 5.2. 



M=2 


a o 


a 1 


a 2 


a 3 


Differential 


. 601 


.425 


.778 


1.20 


Hybrid 


.485 


.485 


.740 


1.05 


M=3 


Differential 


.403 


.285 


.410 


. 525 


Hybrid 


.310 


.310 


.40 


.482 


M=4 


Differential 


.232 


. 164 


.20 


. 231 


Hybrid 


. 172 


. 172 


.20 


.217 


M=5 


Differential 


. 127 


. 090 


. 10 


. 110 


Hybrid 


.092 


.092 


.10 


. 107 



TABLE 5.2 

Error Function Coefficients 
for MF-QAM 
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Theoretical values of SNRyg = E/N 0 = p are given in db in 

-2 “4 *6 

Table 5.3 for BERs of 10 , 10 and 10 due to single bit errors 

for MF16-QAM. (E^^ values are 6 db lower.) 



10~ 2 10~ 4 10 



Differential 17.9 21.9 24.2 

Hybrid 15.9 21.0 23.3 

TABLE 5.3 

Theoretical SNRs in db for BER^ shown 
for Fully Differential and Hybrid MF16-QAM 



b. Double Bit Errors 

Double bit errors are less likely than single bit errors but 
can contribute to the overall BER at the lower values of SNRy B . 

The probability of a double bit error is 

P[B=2) = Pr[ep 2/ e a =0] + Prte^e^l] 

< Pr[ep 2 ] + Pr [e a =l] (5.9) 
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where the last term in the inequality is used as a bound because 
of the probable high correlation between amplitude and phase bit 
errors and is given by (5.8). The double phase error probability 
is determined from (5.5) and (3.35) with (3.27). This is always 
much less than the probability of an amplitude bit error when the 
amplitude spread has been adjusted to a value such that amplitude 
bit errors and single phase bit errors are approximately equal. 
Thus we may approximate the BER due to double bit errors as 



differential coding. 

The effect of double bit errors is to approximately double the 
BER determined by BER,. However, because the Q(x) function falls 
off so rapidly, this can be offset by 0.3, 0.5, and 0.8 db 



Multi-level Amplitude MFOAM Systems 

These systems become quite complex to analyze exactly because 



BER 2 ~ 2Q‘ 1 Pr[e a =l] 



(5.10) 



with 




Q(pa 0 ) for hybrid coding 



(3/2)Q(pa 0 /2*) -(l/2)Q(3pa 0 /2*) for fully 



increases in SNR for BERs of 10 6 ,10’ 4 , and 10’ 6 respectively. 
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of the multiplicity of amplitudes that occur in the complex rotated 
product demodulation of the differential phase terms. However, a 
reasonable hueristic design procedure is to set the differential 
amplitude level such that the distance in the constellation between 
the two smallest amplitudes is equal to the radial distance between 
adjacent phase points on the smallest amplitude ring. This is 
illustrated in Figure 6. The design procedure is to set 



(which leads to a = {1/A <p + (2 N -l)/2 } ') . 

With this constraint, we see that the single bit error probability 
will be bounded by 



which can be found from (4.3) for the amplitude bit errors and from 
(3.35) and (3.27) for the phase bit errors. We now use in addition 
the facts that 




(5.11) 



P[B=1] < Pr [e a =l] + Pr^] 



(5.12) 



AA ' = 2 (2E/N 0 } % f (a , N) 



(5.13) 



with f(a,N) given by (4.11) and A <p = 2tt/ 2 m . Thus, 



Pr[e a =l] = 2Q(AA'/2) 
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and 



Pr[ep, ] = 2Q(M' mjn A0/2) 



with M' mjn = (E mjr /N 0 ) % = A' mir /2 % = AA'/2 V such that 



PrtSpO = 2Q { AA * / ( 2 ( 2^) ) } (5.14) 



It can be noted that the phase errors are 3 db worse than the 
amplitude errors because of the differential coding. However, the 
amplitude errors are the same for all amplitude levels whereas the 

t 

phase errors are much less for the larger amplitude values. Using 
this procedure one obtains the parameters of Table 5.4 for the 
(M,N) hybrids indicated. 



L 




Figure 6 

Ring Spacing to Equalize Amplitude and Phase Bit Errors 
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24.8 


(4,2) 


.248 




. 122 


32.0 


(5,2) 


.152 




.075 


36.3 ,34.3,31 


(5,3) 


. 116 




.058 


38.5 


(6,2) 


.085 




. 046 


40.5 



TABLE 5.4 

SNRs Required for Single BER 
shown for Hybrids from 4 to 8 Bits/Hz. 



The (3,1) 16 -QAM case has been included for comparison. Notice 
how these geometrical arguments lead to results very close to those 
obtained in the previous section by more thourough analysis. (24.8 
vice 23.3 db for 10 6 BER and a of .563 vice .447). This gives us 
considerable confidence in our other results. The (5,2) hybrid 
BERs also seem to be closely verified by our current experiments. 
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VI DISCUSSION AND CONCLUSIONS 



The principal objective of this report was to document the 
theortical predictions for various types of differential coding in 
the frequency domain for Multi-Frequency Modulation (MFM) . In order 
to do this in a systematic fashion, general results were derived 
for probability of symbol error and bit error rate(BER) in Chapter 
II in terms of e, the error pattern vector. The coding we have 
developed partitions a symbol s into phase bits p and amplitude 
bits a. Phase bits are further partitioned into two quad bits g 
that determine the phase quadrant and any remaining phase bits that 
determine tone phase within a quadrant. 

In Chapter III results were presented for coherent phase 
modulation with maximum likelihood decoding and for differential 
phase modulation using the "complex product rotation" algorithim 
followed by maximum likelihood decoding. In order to minimize BER, 
the phase states are Grey coded so that the most probable symbol 
errors result only in single bit errors. 

Differential phase coding and decoding using complex rotated 
product results in a theoretical 3 db loss of performance versus 
SNR (or E^Ng) when considering detection on a tone by tone basis. 
A well known phenomena in differential coding is the increased 
probability of error pairs[8]. This occurs because the noise 
present in a single tone affects the complex rotated product for 
two adjacent symbols. Thus if noise is high in one tone, there is 
an increased probability of two adjacent symbol decoding errors. We 
have not presented a theoretical treatment of this phenomena in 
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this report but will in a subsequent report that treats 
interleaving and forward error control techniques, including 
trellis coding, for MFDM. However, it is mentioned here because in 
addition to the 3 db of performance loss, error pairing is another 
penalty paid for differential coding. 

In Chapter IV results were presented for amplitude coding the 
tones of MFM. Multi-level amplitudes with equally spaced intervals 
were assumed. The performance of maximum likelihood detection of 
the amplitude interval of each tone at the receiver was given. This 
technique is of course incoherent because it does not require 
reference phases for the tones for decoding, but it does require 
knowledge of the absolute gain of the channel at each frequency in 
order to set the decoding thresholds. Consequently, we also 
presented results for differential coding/detection of a single 
amplitude bit. Again, about a 3 db penalty is paid for differential 
coding and the increased probability of error pairing can be 
expected here as well. 

In Chapter V the concepts of phase and amplitude coding were 
combined to create MFQAM constellations. All of these 
constellations are circularily symmetric and because their phase is 
differentially encoded, the decoding is independent of absolute 
tone phase shift thru the channel. Constellations with a single 
differentially encoded amplitude bit and differential phase are 
termed "fully differential" constellations. Constellations with 
absolute level decoding of amplitude bits but differential decoding 
of the phase bits are termed "hybrid" constellations. SNR required 
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to achieve BER of 10 6 were given for a number of the constellations 
in Table 5.4. Constellation design procedures were presented that 
purport to very nearly minimize BER. 

The SNRs (or E^N^ required for differential constellations 
to achieve given levels of BER are about 3 to 5 db higher than is 
required for fully coherent MFM using rectangular (Z 2 ) grid for 
locating the tone real and imaginary parts. The advantage gained by 
differential coding is the removal of any requirement to equalize 
the channel and to determine reference phases for each of the 
tones. In certain applications, the penalities of differential 
coding in SNR will be outweighed by its simplicity of 
implementation, synchronization and robustness in the face of 
changing channel parameters. In other applications, it will not be 
worth the penalty. 

The purpose of this report was to present design procedures 
and performance expectations for frequency domain differentially 
encoded MFM with signalling efficiencies of 2 to 8 bits/Hz. 

One of the newest techniques of error control in single carrier 
digital modulation is trellis coded modulation (TCM) [10] . Our direct 
funded research project is now directed toward applying TCM to MFM 
by trellis coding in the frequency domain. In brief, we believe 
that TCM has the potential for improving our raw Ej/Nj, performance 
by 3 to 4 db. Furthermore, with interleaving it is hoped that the 
effects of differentially induced error pairs can be removed. 

Application of TCM to pure phase coding is believed to be 
relatively straighforward. From an implementaion point of view, the 



50 



MFM receiver must be equipped with a Viterbi decoder. Although 
complex, the procedures are well established to do this and Viterbi 
decoders for TCM have been written for the Motorola 56000 DSP chip 
[11]. Trellis coding the QAM constellations analyzed in Chapter V 
may not be a straightforward application of known results. Previous 
work has concentrated on TCM for Z 2 grids [10] and the conversion to 
our circularly symmetric constellations is not directly apparent. 
This problem is at the top of our current research agenda. 
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